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Demonstration of Low Emittance in the Cornell Energy Recovery Linac Injector 
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Dunham, Francisco Gonzalez, Georg Hoffstaetter, Siddharth Karkare, Hyeri Lee, Heng Li, 
Yuhn Li, Xianghong Liu, Jared Maxson, Christian Nguyen, Karl Smolenski, and Zhi Zhao 

(Dated: April 10, 2013) 

We present a detailed study of the six- dimensional phase space of the electron beam produced by 
the Cornell Energy Recovery Linac Photoinjector, a high-brightness, high repetition rate (1.3 GHz) 
DC photoemission source designed to drive a hard x-ray energy recovery linac (ERL). A complete 
simulation model of the injector has been constructed, verified by measurement, and optimized. 
Both the horizontal and vertical 2D transverse phase spaces, as well as the time-resolved (sliced) 
horizontal phase space, were simulated and directly measured at the end of the injector for 19 pC 
and 77 pC bunches at roughly 8 MeV. These bunch charges were chosen because they correspond to 
25 mA and 100 mA average current if operating at the full 1.3 GHz repetition rate. The resulting 
90% normalized transverse emittances for 19 (77) pC/bunch were 0.23 =b 0.02 (0.51 =b 0.04) /im 
in the horizontal plane, and 0.14 ± 0.01 (0.29 it 0.02) /xm in the vertical plane, respectively. These 
emittances were measured with a corresponding bunch length of 2.1 ±0.1 (3.0±0.2) ps, respectively. 
In each case the rms momentum spread was determined to be on the order of 10~^. Excellent overall 
agreement between measurement and simulation has been demonstrated. Using the emittances and 
bunch length measured at 19 pC/bunch, we estimate the electron beam quality in a 1.3 GHz, 5 GeV 
hard x-ray ERL to be at least a factor of 20 times better than that of existing storage rings when 
the rms energy spread of each device is considered. These results represent a milestone for the field 
of high-brightness, high-current photoinjectors. 



INTRODUCTION 



The desire for light sources with substantially more 
coherence and brightness has fueled significant interest 
in the research and design of ERLs and free electron 
lasers (FELs). The feasibility of ERL technology has 
already been demonstrated at several laboratories, most 
notably Thomas Jefferson National Accelerator Facility 
(TJNAF) [1 , where energy recovery was achieved for 100 
MeV beams with an average current of up to 9 mA. How- 
ever, in order to design and construct a large scale, high 
energy (GeV) ERL x-ray source, significant advancement 
of both superconducting RF (SRF) cavity technology, as 
well as high-brightness, high current sources has been re- 
quired. Over the last several years, Cornell University 
has played a lead role in the development of both areas, 
and has successfully reached several major milestones to- 
wards the realization of a practical ERL x-ray facility [2 . 

To drive this type of machine requires an exceptional 
electron source producing high-brightness bunches at 
high repetition rates. Traditionally it has been thought 
that the best beam quality was obtained using low duty 
factor normal conducting RF (NCRF) gun based pho- 
toemission sources [SHS], as these devices are capable of 
providing high peak cathode fields. These fields are typ- 
ically in the vicinity of 100 MV/m, though the field at 
the cathode during emission is often significantly lower 
since these devices are usually run off-crest HJ [5] . Due to 
the considerable heat load generated in the cavity walls, 
the cw operation of NCRF sources requires substantial 
lowering of the electric gradient, an approach being pur- 



sued at several facilities [3l [6] . Work started at Jefferson 
lab, and later expanded at Cornell University, shows that 
the combination of a high-voltage DC gun followed im- 
mediately by acceleration with superconducting cavities 
yields beams with single bunch quality rivaling that pro- 
duced by RF guns, but at much higher (GHz) repetition 
rates [7, 8 . In addition, DC guns provide an excellent 
vacuum, allowing for a much wider range of cathode ma- 
terials to be used than in NCRF guns. While SRF guns 
show significant promise for producing high-brightness, 
high current beams, the current state of this technology 
is in the development and testing stage, and the achieved 
beam parameters so far are relatively modest [9 . 

Consequently, a photoinjector using a DC gun has been 
designed, built, and commissioned at Cornell University. 
One of the main goals of this project was to demon- 
strate the ability to produce high average current from 
this source. The Cornell injector has made great strides 
toward this end, having recently set a new record for high 
average current from a photoinjector with cathode life- 
time suitable for an operating facility [10 . Another ma- 
jor goal is the demonstration of low emittance at the end 
of the injector's merger section, where the (relatively) low 
energy beam would be injected into the main ERL linac. 
The results in this work demonstrate that it is possible 
to produce and transport beams from a DC source which 
have emittances at the point of injection approaching the 
diffraction limit for hard x-rays, and which have a bunch 
length and an energy spread within the parameter space 
required by the specifications of a full hard x-ray ERL. 

In general, to achieve the maximum brightness in a 



photoinjector, it is crucial to control both the transverse 
and longitudinal space charge forces, as well as the effects 
of time dependent RF focusing [HI fTTHTS] . Effective emit- 
tance compensation is possible when bunches are created 
with a charge distribution that has predominantly linear 
space charge fields [TTJ [T3HT51 , and if done correctly, can 
lead to final emittances approaching the intrinsic emit- 
tance of the photocathode. One fundamental limit to this 
approach occurs when the amount of charge extracted 
from the cathode nears the virtual cathode instability 
limit. A rough calculation shows that the lowest achiev- 
able emittance then becomes proportional to the square 
root of the bunch charge q [15] : 



MTE 



q- 



^cath 



(1) 



Here MTE and £^cath are the mean transverse energy of 
the cathode and the accelerating field at the cathode, 
respectively. Detailed simulations of well optimized DC 
gun photoinjectors support this square root dependence 
on the bunch charge and the cathode's MTE [3|8]. In this 
paper, we show that the final measured emittance also 
scales in accordance with Eq. ([l]). This represents a key 
step in experimentally realizing the maximum brightness 
limit for photoinjectors. 

The outline of this work is structured as follows. First, 
a general description of the Cornell ERL photoinjector is 
given. This includes a description of the beamline lay- 
out, the relevant accelerating and optical elements, and 
the diagnostic systems used to take our emittance data. 
Next, we describe how to model the dynamics in the in- 
jector using the space charge simulation code General 
Particle Tracer (GPT) [16], and give a verification of 
the GPT injector model against linear optics measure- 
ments. After this, a description of the optimization of 
this model and the process for determining our final op- 
tics settings used in the experiment is given. This is 
followed by the main results of this work. These include 
direct measurement and simulation of both the projected 
transverse phase spaces, as well as the time-resolved hor- 
izontal phase space at the end of the injector merger sec- 
tion. Additionally, the energy spread distribution was 
measured using a single dipole magnet in a separate di- 
agnostic beamline section, providing an upper bound on 
the rms energy spread at the end of the merger. 



THE CORNELL ERL INJECTOR 

Construction of the Cornell injector was completed in 
the summer of 2007. Initial beam commissioning exper- 
iments revealed an issue with charging up of the ferrites 
in the higher-order mode dampers in the injector cry- 
omodule. After this problem was successfully addressed 
[17], beam experiments started in earnest in the spring 
of 2010 and have continued to this date [lOl [T8H24] . In 



that time, significant progress towards meeting the tar- 
get goals of the injector project has been made. Table [l] 
shows these specifications. Of particular interest to this 
work are the specifications for the normalized transverse 
emittance and rms bunch length. We demonstrate later 
in this work that these specifications have been met. 



Description and Layout 



TABLE I. List of injector design specifications and target pa- 
rameters. 



Parameter 


Specification 


Beam energy 


5-15 MeV 


Normalized emittance 


Cn < 0.3 fim 


RMS bunch length 


cTt < 3 ps 


Bunch Charge 


77 (19) pC 


Average Current 


100 (25) mA 



The layout of the Cornell ERL injector is shown in 
Fig. [T] The Cornell injector features two laser systems. 
The primary system is a 1.3 GHz laser producing 520 nm, 

1 ps rms pulses with an average power of up to 60 W [25] , 
and is used for high current experiments. For emittance 
measurements with nonzero bunch charge, we exclusively 
use a 50 MHz system, whose individual pulses have com- 
parable pulse energy and duration to the 1.3 GHz laser. 
This laser system allows us to limit the average electron 
beam power hitting our interceptive emittance diagnos- 
tics. After being generated in one of these two lasers, 
the final laser pulse train can be chopped using a Pock- 
ells cell, and shaped using our temporal shaping system 
[To] . This system consists of four rotatable birefringent 
crystals, which are used to divide the primary laser pulse 
into 16 copies, with tunable relative intensities set by 
their rotation angles. These crystals are typically tuned 
to produce a roughly flat intensity profile, with around 8 
ps rms duration. 

The cathode used for this study was a GaAs wafer 
grown using Molecular Beam Epitaxy on a p-doped GaAs 
substrate. The cathode was heat cleaned to 620 °C for 

2 hours and then activated to negative electron affinity 
using Cs and NF3 via the "yo-yo" process. The doping 
density was 5 x 10^^ cm~^. The top 100 nm was left un- 
doped. The resulting cathode had a quantum efficiency 
of 4%, a mean transverse energy of 90 meV, and a sub- 
picosecond response time at 520 nm. 

The high- volt age DC gun used in these measurements 
is the same one used in previous space charge and emit- 
tance studies [lOl fT8H24] . The gun was operated at 350 
kV for all measurements in this work. The beamline sec- 
tion just after the gun, labeled 'Al' in Fig. [T] houses two 
emittance compensation solenoids and a 1.3 GHz normal 
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FIG. 1. Top view of the Cornell ERL injector. 



conducting buncher cavity. These elements were used 
to compensate the initial emittance blow up near the 
cathode, and to compress the bunch longitudinally be- 
fore further acceleration. Immediately after emittance 
compensation, the bunches were accelerated using the 
five superconducting niobium cavities in the SRF cry- 
omodule, labeled 'A2' in Fig. [l] In addition to increas- 
ing the beam energy, and thus partially freezing in the 
emittance, the SRF cavities were also used to perform 
further emittance compensation and longitudinal com- 
pression via time-dependent transverse and longitudinal 
focusing. Each cavity features a symmetric twin input 
coupler design in order to eliminate any time-dependent 
dipole kick [26l [27] and can be operated with a voltage in 
the range of 1 to 3 MV. For a more detailed description 
of the injector cavities see [28] , 

Just after the cryo module, the beam was passed 
through a four-quad telescope, labeled 'A3' in Fig. [T] 
The beam was then directed into one of several diagnos- 
tic beamline sections. The section most relevant to this 
work is the 'Bl' merger shown in detail in Fig. [2] The 
injector merger section is comprised of a conventional 
three-dipole achromat [29H3T] . This design was chosen 
for its simplicity, and due to the limited space available 
for the injector experiment. The trade off for this ap- 
proach is that while this merger set-up closes the single 
particle dispersion, it does not satisfy the second achro- 
mat condition r]'^^ = for the space charge dispersion 
function [29l|30]. Despite this, both our simulations and 
measurements show that this merger design does in fact 
preserve low emittance for our operating parameters. As 
was anticipated in [30] , this was accomplished by finding 
the correct settings for the four quadrupole magnets in 
the A3 straight section. 

The emittance measurement system (EMS) used for 
projected and time-resolved phase space measurements 
is a two-slit system with no moving parts ^21,. Fig. |2] 
shows the layout of this diagnostic system. In front of 
each 20 /im slit is a scanner magnet. Each scanner mag- 
net consists of a pair of air core correcting coils with equal 
and opposite field polarity and negligible sextupole field 
component. The resulting effect of the scanner magnet 
is to translate the beam transversely without imparting 



any angle to it. In practice the coil pairs in each scan- 
ner magnet cancel each other to better than a few per- 
cent [21 . For projected phase space measurements, the 
beamlet coming through both slits was collected using 
the Faraday cup at the end of the merger section. For 
time-resolved horizontal phase space measurements, the 
beamlet was passed through a horizontal deflecting cav- 
ity [32 in order to resolve the time axis of the beam on 
the viewscreen at the end of the merger section [24] . For 
a more detailed description of the EMS, refer to [2TJ[24]. 



The GPT Injector Model 

The 3D space charge code GPT was used extensively 
in this work. To model space charge effects, GPT uti- 
lizes a 3D non-equidistant mesh solver [33l |34]. Addi- 
tionally, GPT allows users to define their own custom 
optical elements, as well as position and overlay electro- 
magnetic field maps in 3D space. These features pro- 
vided sufficient versatility to accurately model our ma- 
chine, where the fields of several elements overlap. All 
of the beam line elements relevant for the space charge 
simulations in this work have been modeled using real- 
istic field maps. POISSON-SUPERFISH [35] was used 
to generate 2D cylindrically symmetric fields specifying 
Erir^z) and Ez{r^z)^ as well as Brir^z) and Bz{r^z), 
for the high- volt age DC gun and emittance compensa- 
tion solenoids respectively. The on-axis fields for these 
elements are shown in Fig. 3(a)[ and |3(b)[ 

In order to efficiently and accurately describe the in- 
jector dipoles and quadrupoles, we created custom GPT 
elements which generate 3D fields using an off-axis field 
expansion of ID field data. To create the ID dipole and 
quadrupole field data, the full 3D fields for each type of 
element were computed in Opera-3D [36 . From these 
fields the quantities By{r = 0, z) and dByir — {)^z)/dx 
were extracted from the dipole and quadrupole fields re- 
spectively. Fig. |3(c)| and |3(d)| show the ID field data 
used for the dipoles and quadrupoles in the injector. Our 
custom GPT rectangular dipole element uses an off-axis 
expansion of the fields given by [37, : 



B, ~ 0(4) 
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FIG. 2. Top view of the Bl injector merger section showing the emittance measurement system. 



By 



Bz = y 



Bo{z) 
dBo 



dz 



0{4), 



0(4), 



(2) 



to model the the higher order components of the dipole 
field. In this expression Bq = By{r = 0,z). This ex- 
pansion assumes that the particles do not see the fringe 
fields on the lateral sides of the magnet. This is true 
for the particle trajectories and magnets in the injector, 
where the maximum simulated rms beam size through 
the dipoles was < 3 mm (see Fig. 13(b)), and the dipole 
width was 25 cm. Similarly, the fields for the quadrupoles 
were computed with an off-axis field expansion [37]: 
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0(5), 
■0(5), 



B, = xyG{z) + 0(4). 



(3) 



Here the term G{z) = dBy/dx{r = 0,z). To verify 
Eqs. ([2]) and ([s]), single particle tracking through the 
fields created by our custom elements was compared to 
tracking using the full 3D field maps. Excellent agree- 
ment was found in both cases. Additionally, the custom 
elements proved significantly faster because they do not 
require look-up of 3D field arrays. 

All RF cavity fields were generated using the eigen- 
mode 3D field solver in GST Microwave Studio (MWS) 
[38 . The buncher cavity was modeled using a 2D cylin- 
drically symmetric map specifying: £^^(r, 2:), Ez{r^z)^ 
and Ho{r^z). The on-axis field map for the buncher is 
shown in Fig. 3(e)[ 

Previous work demonstrates that asymmetric focusing 
of the bunch near the input power couplers of the accel- 
erating cavities is significant and can lead to asymmetric 



horizontal and vertical emittances [2T. To address this 
issue, we generated full 3D field maps for the accelerating 
cavities which incorporate the beam running conditions 
following our method outlined in [39]. Fig. El shows the 
3D cavity model used in MWS for the accelerating cavi- 
ties. The procedure for correctly constructing the fields 
in the coupler and cavity requires two sets of MWS so- 
lutions. Each set of fields was created by terminating 
the input coupler line in the MWS model with either an 
electric or magnetic wall boundary condition [39 . From 
these solutions, traveling waves carrying power into and 
out of the cavity through the couplers were constructed, 
scaled, and shifted in phase to match the actual run- 
ning conditions in the injector. In order to further limit 
the beam power deposited in our interceptive EMS, the 
pulse train from the 50 MHz laser was chopped using a 
Pockells cell. The resulting beams typically had currents 
on the order of a micro-amp or less. In generating the 
field maps for the accelerating cavities, this amounts to 
effectively having zero current. The parameter which de- 
termines how the cavity fields depend on the operating 
parameters is the reflection coefficient [39] : 

fef + ^(i?/g)QLe-*^° +ztanV^^ 



1 + ff (i?/g)QLe-^^° + ^ tan V^^ 

Here I^ and 0o are the average beam current and phase 
of the beam with respect to the cavity fields. The rest 
of the parameters in this expression describe the proper- 
ties of the cavities: P is the coupling parameter, Q is the 
intrinsic cavity quality factor, R is the shunt resistance, 
Vc is the cavity voltage, Ql = Qo/{l -\- P) is the loaded 
quality factor, and tan?/;' is the cavity detuning param- 
eter. In the zero current limit (with the cavity tuned to 
resonance), the reflection coefficient reduces to 
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(a) On-axis electric field in the high voltage DC gun. 




-0.1 0.1 0.2 

z(m) 

(b) On-axis magnetic field for the Al solenoids. 
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(e) On-axis electric field for the buncher. 
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FIG. 3. On-axis electric and magnetic fields for: (a) the high voltage DC gun at 350 kV, (b) emittance compensation solenoid 
at 3.75 A, (c) the A3 and Bl merger dipoles, (d) the A3 and Bl merger quads, (e) the buncher cavity at 60 kV, (f) the SRF 
cavity at 1 MV. 



This implies that the fields in the cavity and coupler coax 
depend only on the amount of coupling. The coupling 
factor /3 is determined by how far the inner coupler an- 
tennae are retracted from being fiush with the beam pipe. 
For our emittances measurements, the couplers were fully 
retracted (zero current setting). By fully retracting the 
coupler antennae in the MWS model, and generating two 
set of solutions for both boundary conditions on the end 
of the coupler coax, we created one set of complex 3D 
electric and magnetic field maps for the SRF cavities. 
Fig. |3(f)| shows the resulting on-axis electric field. We 



point out that simulations subsequently showed that the 
asymmetric emittances caused by the RF quad effect in 
the cavities could be successfully remedied by appropri- 
ate choice of magnetic quadrupole focusing downstream. 

After completing the GPT physics model of the injec- 
tor, we developed a user interface between the real ma- 
chine and its GPT counterpart. Named the "Virtual Ac- 
celerator GUI" , this program was designed to provide a 
single interface between the corresponding optics settings 
in the EPICS control system of the real machine, and the 
stand alone GPT code. Additional features include the 




(a) Exterior view of the MWS injector cavity model, (b) Cutaway view of the MWS injector cavity model. 



FIG. 4. (a) The Microwave Studio model of ERL injector cavity: (a) the cavity and coupler exterior, (b) cutaway view of the 
same model showing the inner coupler antennae. 



ability to save and load optics settings and simulation 
results to and from file, the ability to load injector set- 
tings from the machine and independently adjust them 
in simulation, as well as the ability to visualize all rele- 
vant simulation data. A screen shot of this application is 
shown in Fig. [5] In constructing this program, a master 
GPT input file was created which included not only the 
optical elements described in this section, but also simu- 
lation output screens at all of the corresponding locations 
of the BPMs, viewscreens, and emittance measurement 
systems in the injector. The result was a nearly one-to- 
one simulation counterpart to the real machine. With 
this, we were able to use GPT in a more useful and real- 
istic way, with simulations often guiding experiments in 
near real-time in the control room. 



MEASUREMENTS 

All of the measurements in this work fall into one of two 
categories: measurements performed at near-zero bunch 
charge for verification and calibration purposes; or phase 
space measurements of space charge dominated bunches. 



Measurements at Near-Zero Bunch Charge 

The measurements presented in this section include: 
comparison of difference orbits (linear optics) in the in- 
jector with the GPT model including the effects of the RF 
input couplers, measurement of the beam size envelope 
along the injector and its verification with simulation, 
and calibration of the EMS and analysis procedures by 
comparing the emitance computed from the direct mea- 
surement of the projected transverse phase spaces in the 
merger and the emittance measured using a solenoid scan 
in the gun vicinity. In order to accurately perform differ- 



ence orbit measurements, the BPM system needed to be 
corrected for its non-linear response, the procedure for 
which is presented below. 



BPM Correction Procedure 

The injector BPMs consist of four striplines, as seen in 
Fig |6(a)[ To model the system we make two assumptions: 
(i) both the beam pipe and striplines are assumed to be 
infinitely long perfect conductors connected to ground; 
and (ii) the beam is taken to be an infinite line charge at 
the position r^ = (x^, 7/5). The first assumption implies 
that the potential must vanish at the beam pipe. This 
is accomplished by placing an image line charge with op- 
posite charge density at r 
electric field everywhere is: 



{R^/rl)rij. The resulting 



E(r,r5) 



A 



27reo 



r - r5 r - {R'^/rl)n 



n 



{Ryr 



)n 



From the electric field, the surface charge density on the 
beam pipe and striplines can be computed using a = 
— eoE(R, r^) -n, where n is the normal vector to the beam 
pipe surface. Since this model assumes the beam pipe 
is a perfect conductor, the field is perpendicular to the 
surface so that E(R, r5) • n = |E(R, r5)|. In cylindrical 
polar coordinates, the surface charge density takes the 
form 



cT{R,0,rb,Ob) 



Jrl/R') - 1 
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The angles 6 and O^ are defined in Fig 6(a) The signal 
i^^ stripline is defined as the fraction of the 
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surface charge density found on that stripline: 
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Beamline and Save/Restore- 



Select Beamline Tvpe_ 

ala2a3bl 

II Load Saved Settings 
I Save Settings 

I Load Machine State 
earn Element Setting 



-GPT Simulation Settings 



nps 


20000 [#] 


Qtot 


-??[pC] 


space.charge 


1[0/1] 


xyrms 


0,4700 [mm] 


trms 7,9610 [ps] 


xoffset 


0[mm] 


yoffset 


0[mm] 


giobaLphase 


[deg] 


dtmax 


100 [ps] 


FinalCamma 


^^^[N/A] 


couplers 


1[0/1] 




Run GPT Simulation- 



Pliase GPT 



<- or -> Just Set Energy 



- 1 particle, no space cinarge 

- Accurate witiiin few tenths of degree 

- Phases a cavity on crest, then sets 
it off crest (as desired) before 
phasing the next cavity 



- no space charge 

- Assumes user l^nows what he is 
doing and just runs the code 
through A2 to set the final energy. 






Element 


1 PV Name 




\ I Position 1 




1 


GAIGVHOI 


Gun Voltage 


349.3000 


349.3000 (k.V) 


349.3000 (kV) 





2 


MAlSLAOl 


Current 


-3.7000 


-3.7069(A) 


0.0309 (T) 


0,3030 


3 


RAlCTBOl 


Cavity Voltage 


60 


60 (kV) 


0.9569 (MV/m) 


0,7140 


4 


RAlCTBOl 


On-crest Cavity Phase 


175.8100 


175.8100 (deg 


148.4674 


0,7140 


5 


RAlCTBOl 


Relative Cavity Phase 


-90.0000 


-90.0000 (deg 


-90.0000 


0,7140 


6 


MA1SLA02 


Current 


2.3000 


2.3023(A) 


-0.0192 (T) 


1,1280 


7 


RA2CTC01 


Cavity Voltage 


1491 


1491(kV) 


12.7983 (MV/m) 


2,0470 


8 


RA2CTC01 


On-crest Cavity Phase 


-61.2000 


-61.2000 (deg 


9.1949 


2,0470 


9 


RA2CTC01 


Relative Cavity Phase 


-10.0000 


-10.0000(deg 


-10.0000 


2,0470 


10 


RA2CTC02 


Cavity Voltage 


1900 


1900 (kV) 


16.3090 (MV/m) 


2,8330 


11 


RA2CTC02 


On-crest Cavity Phase 


43.2000 


43.2000 (deg 


356.1229 


2,8330 


12 


RA2CTC02 


Relative Cavity Phase 


-7 


-7.0000 (deg 


-7 2,8330 


13 


RA2CTC03 


Cavity Voltage 


1.3860e+03 


l,3S60e+03(kV) 


11.8970 (MV/m) 


3,6950 


14 


RA2CTC03 


On-crest Cavity Phase 


161.0900 


161.0900 (deg 


256.8649 


3,6960 


IS 


RA2CTC03 


Relative Cavity Phase 


-3,5621e-06 


-3.6621e-06(deg 


-3,5521e-06 


3,6950 


16 


RA2CTC04 


Cavity Voltage 


1386 


1386 (kV) 


11.8970 (MV/m) 


4,4820 


17 


RA2CTC04 


On-crest Cavity Phase 


97.1322 


97. 1322 (deg 


284.5580 


4,4820 


18 


RA2CTC04 


Relative Cavity Phase 
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FIG. 5. Screenshot of the Virtual Accelerator GUI. 



Here the angle Os is the angle subtended by each stripline. 
Performing the integration yields: 
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ei+9s/2 



, (7) 



./2 



where Oi G {0, 7r/2, tt, 37r/2}. 

In order to invert the BPM signals and obtain the beam 
position, the signals from Eq. (ItI are fit to the injector 
BPM signals using a x^-minimization with the beam po- 
sition as the fit parameters. To verify this procedure, a 
pair of upstream horizontal and vertical corrector mag- 
nets was scanned in a grid pattern and the response on a 
test BPM was measured. Fig |6(b)| shows the comparison 
of the standard linear BPM model (blue), and non-linear 
model given by Eq. ^ (red). The inclusion of this model 
effectively extended the workable range of the BPMs in 
the injector by roughly a factor of two. This increased 
range made the use of the BPMs in response measure- 
ments significantly more robust. 



Difference Orbits and Coupler Effects 

To verify each injector beamline element and its cor- 
responding GPT model, linear optics response measure- 
ments have been performed. The transverse dynamics 
were verified by changing the initial position of the beam 
on the cathode or kicking the beam with a corrector mag- 
net and recording the change in position on all down- 
stream BPMs. This was repeated for each type of ele- 
ment in the injector, starting with the gun and moving 
downstream turning on elements one by one and com- 
paring the resulting response function to GPT simula- 
tions. Fig. 7 (a) I shows an example response measurement 
and corresponding GPT comparison. For this measure- 
ment, the first pair of horizontal and vertical correctors 
in the Al section were scanned and the response through 
the straight portion of the injector recorded (with all 
quadrupoles off). Time of flight difference orbits were 
also measured by adjusting the laser phase ±60 degrees 
relative to the cavity phases. Excellent agreement with 
the GPT model was obtained using all BPMs, including 
those in the merger. 

Difference orbits were also used to verify the 3D RF 
field maps used to model the cavities and fields near the 
input power couplers. Simulations show that asymmet- 
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(a) Description of non-linear BPM model parameters. 
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FIG. 6. Non-linear BPM model description and verification: (a) shows the relevant parameters for the model, while (b) shows 
the comparison of the standard linear BPM position calculation (blue) and positions computed with a non-linear correction 
(red) from a square grid scan of an upstream horizontal and vertical corrector pair. 
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(b) Asymmetric response through the second cavity due to coupler 
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FIG. 7. Response measurements: (a) the response from the set of correctors through first SRF cavity. The cavity was set to 1 
MV on-crest. (b) The response asymmetry due the coupler fields in the second cavity as a function of cavity phase. The cavity 
voltage was 1.5 MV. The dashed line shows the expected response from a cylindrically symmetric, or ID field map model of 
the cavity. 



ric focusing from the couplers is worse when a low energy 
beam passes through the coupler fields before being ac- 
celerated [39 . Thus to more clearly measure the effects 
of the couplers, we turned off all of the SRF cavities ex- 
cept the second one, which has couplers at the entrance of 
the cavity, as seen by the beam. A square grid of angles 
was scanned using the last pair of horizontal and vertical 
correctors just before the entrance to the cryomodule, 
and the resulting response pattern was measured on a 
downstream BPM. This was repeated at multiple cavity 



phases shifted relative to the on-crest phase. By taking 
the ratio of the change in position in x to the change in 
y^ the asymmetry in the response though the cavity was 
computed. Fig. |7(b)| shows the comparison of the x to y 
response aspect ratio measured in the injector and com- 
puted in GPT. The agreement is quite good except for 
the point where the response in both planes goes through 
zero. With these measurements, we are confident in our 
ability to include the 3D focusing effects of the cavity 
input couplers. 



Alignment 

Previous work has shown j8l [2TJ [24] that good ahgn- 
ment through each optical element is required to dimin- 
ish emittance growth, and indeed alignment of the beam 
through the gun, emittance compensation section, and 
SRF cavities proved very important for obtaining the 
low emittance results presented here. In order to arrive 
at these results, a methodical element by element align- 
ment procedure was developed. The benefit of such an 
approach was that after a thorough execution of the fol- 
lowing procedure, additional alignment work was kept to 
a minimum on subsequent experiments. The first step 
of this procedure was to center the laser spot on the 
cathode. To do so, the spot was scanned both horizon- 
tally and vertically to form a grid of positions. With 
the first solenoid off, the corresponding beam position 
was recorded on the viewscreen after the gun, see Fig. |8] 
Fitting the response data to an off-axis expansion of the 
gun focusing allowed us to the determine electrostatic 
center of the gun/cathode to within 50 fim. It should be 
noted that in order to achieve a good cathode lifetime, 
off-center laser spot operation is required to minimize ion 
back-bombardment [10 . However, we found that offset- 
ting the laser spot by 3 mm and using a corrector pair 
to bring the beam back through the center of the 1st 
solenoid did not degrade the beam emittance by more 
than 5%. 

After aligning to the gun, the beam was then aligned 
in the buncher cavity. To do so, the gun was set to 350 
kV, and the first and second solenoids were degaussed 
and turned off. The buncher cavity was turned on at 
50 kV and the two energy zero-crossing phases deter- 
mined. In order to keep the transverse beam size small, 
the cavity phase was set to the de-bunching zero-crossing 
value, in order to provide focusing from the buncher. The 
use of the zero-crossing phase also eliminated the effect 
of dispersion due to the combination of unwanted stray 
fields and low beam energy. The beam position on the 
second viewscreen was recorded with the cavity turned 
off and then turned on. The initial position offset going 
into the cavity field region was then found by fitting the 
beam transfer matrix from the corrector coil pair just be- 
fore the buncher to the viewscreen after the cavity. The 
transfer matrix was computed from the on-axis electric 
field map shown in Fig. 3(e) [ using the method derived in 
[39] . The position offset in the buncher was then com- 
pensated by adjusting the corrector coils just before it. 
Using this technique, we were routinely able to align the 
beam through the center of the buncher to within 20 /im. 

Next, the orbit was aligned through the first two SRF 
cavities. Each cavity was separately turned on to 50 kV 
and set to the de-bunching zero-crossing phase just as 
with the buncher. Once the correct phases were found, 
the beam position on the A3 viewscreen was recorded for 



three different settings: both cavities off, and then each 
cavity on separately. After recording the beam position 
on the A3 viewscreen for each setting, the response func- 
tions from the last two pairs of horizontal and vertical 
correctors before the cryomodule were measured. From 
this set of response measurements, the corrector settings 
were determined that would place the beam at the same 
spot on the A3 viewscreen for all three cavity settings. 
This process produced an orbit which did not change 
position on the A3 viewscreen to within roughly 50 jam 
when the first two cavities were toggled on and off. 

Finally, the solenoids were aligned. The alignment of 
the buncher and first two SRF cavities fixed the settings 
of all the available corrector coils in the Al section. Con- 
sequently, the solenoids had to be physically moved to 
align their magnet centers with beam orbit. For the 
solenoids, both their offset and angle in the horizontal 
and vertical planes were found by performing a current 
scan of each magnet, recording the response on a down- 
stream viewscreen, and fitting the data using the transfer 
matrix of the solenoid [23l[40]. The physical adjustment 
of the solenoid positions and angles was greatly aided by 
the incorporation of alignment motors in the design of 
the solenoid magnet support structure. At the comple- 
tion of the final alignment measurements, the transverse 
offsets of the solenoids were aligned to within roughly 50 
/im, and the transverse angles to within 0.2 mrad. 

Alignment of the orbit through the optical elements 
in the A3 straight and Bl merger section was achieved 
by flattening the BPM readings in these sections. To 
check the overall alignment once the orbits for emittance 
measurements were set up, a special laser mask with a 
regular grid of 100 /im holes spaced 0.75 mm apart, was 
placed in the laser path. Fig. [9| shows the initial grid 
pattern and the measured grid pattern in the Bl section. 
To generate this image the buncher was purposefully set 
to give a longer bunch length in order to exaggerate the 
time dependent RF focusing from the SRF cavities. In 
this image, the center spot being circular and the other 
spots pointing towards the center indicates that the beam 
is aligned reasonably throughout the injector. Also, the 
lack of curvature to each of the spots/lines indicates that 
aberration effects are minimal. 



Beam Sizes and Thermal Emittance 

Before measuring emittance with nonzero bunch 
charge, we calibrated our emittance measurement sys- 
tem and emittance analysis scripts by measuring the 
emittance at near-zero bunch charge {q < 0.03 pC). A 
baseline thermal emittance was measured after the gun 
and before the cryomodule by scanning the current of 
the first solenoid and measuring the beam spot size on a 
viewsreen downstream. By computing the linear trans- 
fer matrix through the combined gun and solenoid fields. 
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FIG. 8. Layout of the Al emittance compensation section. 




(a) Initial grid of laser spots as measured (b) Resulting beam image on the 

on the laser ccd. viewscreen at the end of the Bl merger. 



FIG. 9. Alignment check using a grid of laser spots (a) and the resulting beam image on the viewscreen at the end of the Bl 
merger section (b). 



the emittance and initial rms beam spot size were found 
using the method in [23]. Fig. |lQ(a) and |lQ(b)| show 
the solenoid scan data and fitted curve for the cathode 
used in this work. The resulting horizontal and vertical 
emittances measured with the solenoid current scan were 
0.12 ±0.01 /im and 0.11 ±0.01 /im respectively. To check 
the calibration of the EMS in the merger section, the pro- 
jected horizontal and vertical phase space, as well as the 
horizontal time-resolved phase space were measured. For 
these measurements, the 19 pC/bunch injector optics set- 
tings were used (see Table [ll|, however the bunch charge 
was reduced so that space charge effects were negligible. 



Fig. 10(c)| and Fig. |10(d)| show the projected emittance 
measured in the Bl merger section with a beam momen- 
tum of roughly 8 MeV/c. The estimated systematic error 
in the calibration of the merger EMS system was less than 
7%. The horizontal and vertical emittances from these 
measurements were 0.11 ± 0.01 and 0.12 ± 0.01, which 
agree with the solenoid scan results to within the esti- 
mated error in both measurements. The same value for 
the horizontal projected emittance, 0.11 ± 0.01 /im, was 
measured in the merger section using the time-resolved 
EMS. These measurements not only verified the EMS di- 
agnostics and analysis procedures, but also provided an 
additional check of the orbit alignment. 
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(c) Measured horizontal phase space in the merger section. 
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(d) Measured vertical phase space in the merger section. 



FIG. 10. Projected emittance measurement at the cathode using a solenoid scan (a-b), and corresponding measurements in 
the merger section (c-d). Both the horizontal and vertical emittance measured in the merger section agreed to within 9% of 
values measured at the cathode. The colormap and normalization in (c-d) is used for all subsequent phase-space plots in this 
work. The estimated error in for these emittance values was ±0.01 /xm. 



As a final check of the optics settings in the machine 
and simulations, we measured the transverse rms beam 
sizes at several locations along the injector with near-zero 



bunch charge. Fig. 11 shows the comparison of the sim- 
ulated and measured rms spot sizes. The optics settings 
were the same as those used in the EMS calibration mea- 
surements. The measured values were computed from 
images of the beam on the Al and A3 viewscreens, and 
from the phase spaces measured in the merger section 
shown in Fig. 10(c) and Fig. 10(d) The systematic un- 



certainty in these measurements due to the viewscreen 
calibration and set-up resolution was estimated to be 
less than 5% for the direct viewscreen measurements. As 



Fig. 11 shows, excellent agreement between GPT and the 
measured beam sizes was found. 



Measurements with Space Charge 

Two main data sets were produced for this work: one 
at 19 pC per bunch, and one at 77 pC. These correspond 
to 25 mA and 100 mA average current when operating 
at the full 1.3 GHz repetition rate. Each data set con- 
sists of a measurement of the projected horizontal and 
vertical phase spaces, the time-resolved horizontal phase 
space, and the energy spread distribution. All data was 



12 






-^x 


GPT 




a 

■ V 


GPT 


o 


a 

X 


data 


D 


^ 


data 




14 



s(m) 



FIG. 11. Comparison of simulated and measured rms spot sizes along the beam line. 



taken at the end of the merger section except the energy 
spread data, which was measured using the A4 straight 
section and C2 bend section. From the projected phase 
space the horizontal and vertical emittance as a func- 
tion of beam fraction was computed. Similarly, from the 
time-resolved phase spasce data, the slice emittance was 
computed as a function of beam fraction, as well as the 
current profile along the bunch. Refer to Appendix A 
for the emittance definitions used to characterize non- 
Gaussian phase spaces. 



Injector Settings and Simulation Parameters 

To arrive at the final optics used for these experiments, 
optimizations of the GPT model were carried out using a 
multi-objective genetic algorithm [71 [8]. In general, each 
optimization was run with two competing objectives (e.g. 
minimizing the emittance at the location of the merger 
EMS and maximizing the bunch charge), while varying 
the optics settings (e.g. solenoid, rf, and quad settings). 
Upon convergence of the optimizer, this produced an op- 
timal front for the two objective variables. A complete 
list of the parameters varied in the optimizer can be found 
in the first and second columns of Table [ill Note that for 
all optimizations, the gun voltage was fixed at 350 kV, 
and the beam energy was constrained to be < 8 MeV 
to reduce neutron production from the tungsten slits in 
the EMS. The simulated temporal laser distribution was 
fixed to be roughly a flat-top with 8 ps rms length, and 
was generated by adding 16 Gaussian pulses in accor- 
dance with the temporal laser shaping system used for 
the injector [25]. The transverse laser profile was a Gaus- 
sian truncated at 50% intensity with the resulting rms 
size varied in the optimizations. 

Each simulation was run with a set of constraints which 
ensured the physicality of the results and pushed the op- 
timizer to explore regions of the variable space relevant to 



achieving the injector design goals. The two most impor- 
tant constraints were the rms bunch length: at ^ 2 ps, 
and the rms energy spread as ~ 0.1 to 0.2 %. Note that 
these values are slightly more stringent than the quoted 
design goals in Table [l[ 

Optics solutions from the last set of optimizations were 
loaded into the GPT virtual accelerator GUI, and then 
tested in the injector. Doing so led to the recognition 
of a common feature to all of the optimized solutions: 
the optimizer always focused the beam through a waist 
at exactly the position of the EMS in the merger. As 
the experiment proceeded, another general trend was ob- 
served using the virtual accelerator GUI. Optimized so- 
lutions which kept the beam sizes small, particularly in 
the straight section and merger, gave better measured 
emittance results. The settings used in the final mea- 
surements presented here were both derived from one 
optimization solution for 50 pC bunch charge, as this op- 
tics setting kept the beam sizes reasonably small through 
the entire injector. Using this parameter set as a start- 
ing point in our virtual accelerator interface, the bunch 
charge was reduced from 50 to 19 pC while adjusting the 
magnet and buncher settings to compensate for the re- 
duced space charge effects, as well as scaling the laser 
spot diameter in accordance with ax^y ex y^. This al- 
lowed us to keep the simulated spot sizes small through 
the injector, while also maintaining the location of the 
beam focus at the position of the merger EMS. The pro- 
cedure was then started over, raising the charge from 19 
to 77 pC. This time the phases of the first two SRF cav- 
ities, and the voltage and phase of the last SRF cavity 
were also adjusted in order to maintain small emittance 
values at the merger EMS. 

These settings were then loaded into the injector and 
the measured projected emittance was minimized by 
scanning both solenoid currents and adjusting the inten- 
sity cut off value in the measured transverse laser profile. 



Fig. 12 shows the measured laser profiles for used in the 
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final measurements and the corresponding profiles used 
in the final GPT simulations. The final solenoid currents 
used in the injector were within 3% of the simulation 
values. The quads in the Bl section were also adjusted 
slightly for both optics settings, but kept within 4% of 
the simulations. These slight adjustments to the simu- 
lated injector settings are believed to be a consequence 
of hysteresis effects in the magnets, as well as error in 
the calibration factors used to convert machine param- 



eters to simulation parameters. Fig. 13 shows the rms 
beam sizes, projected horizontal and vertical emittances, 
kinetic energy, and bunch lengths computed using the fi- 
nal simulation optics values. Table [ll] shows injector set- 
tings and parameters used in measurements. The beam 
kinetic energy measured after the cryomodule was 7.5 
and 7.7 MeV for the two bunch charges respectively. 



Projected Emittance Results 

As previously discussed, measuring low projected emit- 
tance after the merger section that scales according 
Eq. ([l]), and meets the design specification of the injec- 
tor was one the main goals of this work. Tables |III(a)| 
and |III(b")] show the best projected emittance data from 
measurement as well as the corresponding GPT simula- 
tion values. The measured emittance data was processed 
with removal of a near constant background via an au- 
tomatic bias determination routine similar to the meth- 
ods described in [HI [41] . The processed data was then 
used to determine the 100% beam emittance, as well as 
to generate the emittance vs. fraction curve, defined in 



Eqs. ( 10 ) and ( 11 ), and the corresponding core emittance 



and core fraction, defined in Eq. (12). These curves are 



shown for the horizontal and vertical projected phases at 
19 (77) pC/bunch in Fig. 14 All of these procedures were 



automated and available to operators in the control room 
after each emittance measurement scan (lasting typically 
several seconds). Further details of the data processing 
and experimental procedures can be found in [42 . The 
measured 19 (77) pC /bunch horizontal and vertical pro- 
jected 100% emittances agreed with the GPT model to 
within 6 (5) % and 25 (8) %, respectively. Similarly, the 
measured horizontal and vertical 90% emittances agreed 
with GPT to within 21 (16) % and 27 (16) %, respec- 
tively. We point out that the measured horizontal and 
vertical 100%, 90%, and core emittances obey the ex- 
pected scaling law e^ ex ^yq. Also of note is the fact the 
horizontal core emittance for 77 pC meets the injector 
design specification for an ERL. In the vertical plane, 
both the 90% and core emittance meet this specifica- 



tion. For comparison purposes. Fig. 15(a) and Fig. 15(b) 



show the measured and simulated phase spaces after the 
merger for both the horizontal and vertical planes with 
near-zero, 19, and 77 pC/bunch. 



Time-resolved Phase Space and Energy Spread Results 

In order to satisfy the injector design requirements, it 
was important to verify that the emittance values were 
measured with an acceptable bunch length {at < 3 ps). 
The rms bunch length was computed from the instan- 
taneous current of each bunch measured with the time- 
resolved merger EMS. Fig. |16(a) and |16(b)] shows both 
the measured and simulated bunch current for the 19 
pC/bunch and 77 pC/bunch data, respectively. The rms 
bunch lengths for the 19 (77) pC per bunch settings were 
measured to be 2.1 ±0.1 (3.0 ±0.2) ps, respectively, while 
GPT gave bunch lengths of 2.2 (3.1) ps, respectively. The 
agreement between measurement and GPT was within 
5% in both cases. As Fig. |16(a)| shows, the qualitative 
agreement between data and simulation was good for the 
19 pC /bunch measurement. The difference in the over- 
all scaling between the measured and simulated data for 
this setting is due to the normalization of the data to the 
bunch charge. In the 77 pC /bunch case, the qualitative 
agreement between measurement and simulation, shown 
in Fig. 16(b), was excellent. 

In addition to measuring the bunch length, the time- 
resolved emittance was measured both because it is of 
interest to FEL applications and to elucidate the char- 



acter of the emittance growth in the merger. Fig. 17 
shows the core and 90% emittance for both the 19 and 77 
pC measurements. For the 19 pC data, both emittances 
are relatively constant over the bunch length. Similarly, 
for the 77 pC data, the core emittance is constant over 
the majority of the bunch length. Also important is the 
fact the core emittance for this data is below the design 
specification for the injector. The time-resolved emit- 
tance measurements provide a very elegant way of view- 



ing the resulting phase space distributions. Fig. 18 shows 
the three-dimensional representation of the time-resolved 
phase space for both data sets. The 3D representation 
demonstrates that the z-shaped features seen in the pro- 
jected emittance in Fig. |15(a) are actually a real effect 
formed along the time axis. 

The last quantity measured was the rms energy spread. 
To do so, the beam was sent through the A4 straight sec- 
tion, followed by a single dipole and viewscreen in the 
C2 section (see Fig. IT]). Before entering the dipole, the 
beam was clipped by passing it through a crossed pair 
of emittance measurement slits. The emittance measure- 
ment scanner magnets in this section were set so that 
crossed slits selected out a beamlet from the centroid of 



the horizontal and vertical phase spaces. Table IV shows 
the simulated and measured rms energy spread in the 
straight section, as well as simulated values in the Bl 
merger. The measured values were computed from the 
ID dimensional energy spread distribution obtained on 
the viewscreen in the C2 section using a 10% thresh- 
old to remove background noise. The values are slightly 
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(a) Laser spot for the 19 pC/bunch data set. 
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(b) Laser spot for the 77 pC/bunch data set. 





(c) Radial laser intensity for the 19 pC/bunch settings. 
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(d) Radial laser intensity for the 77 pC/bunch settings. 



FIG. 12. Verification of the initial transverse laser spot in the injector and GPT input file: (a) and (b) show the measured 
laser spot on a ccd camera as they would appear on the cathode. The uncertainty in the rms spot sizes is (a) ±0.01 mm, and 
(b) ±0.02 mm. The plots in (c) and (d) show the corresponding measured radial laser intensity (blue) and matching truncated 
Gaussian (red) used in simulation. 



smaller than simulation, which is likely due to the fact 
that we were only measuring the energy spread of a sin- 
gle transverse beamlet, while the simulated values are 
computed from the entire beam distribution. While we 
did not measure the energy spread directly in the merger 
section, the agreement found between measurement and 
simulation for emittance and bunch length lead us to con- 
clude that the values measured in the straight section at 
least provide an upper bound on the energy spread in the 
merger, following the same trend found in the simulation 
data. 



CONCLUSION AND DISCUSSIONS 

A comprehensive model of the Cornell ERL injec- 
tor has been constructed using the space charge code 
GPT. After verifying the accuracy of the GPT model 
against linear optics measurements in the injector, multi- 
objective optimizations of the model were carried out in 
order to find optics settings with which to measure low 
emittance after the merger section in the injector. In ad- 
dition, a user interface between the GPT code, the opti- 
mizer solutions, and the injector was developed. This in- 
terface provided visualization of relevant simulation data 
in one-to-one correspondence with measured data, and 
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(a) RMS beam sizes for 19 pC/bunch. 
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(b) RMS beam sizes for 77 pC/bunch. 
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(c) Projected emittance for 19 pC/bunch. 
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(d) Projected emittance for 77 pC/bunch.. 




2 4 6 8 10 12 

s(m) 

(e) Kinetic energy and bunch length for 19 pC/bunch. 
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(f) Kinetic energy and bunch length for 77 pC/bunch. 



FIG. 13. Simulation data for the 19 pC/bunch (left) and 77 pC/bunch (right) injector settings: (a-b) show the rms beam size 
along the injector, (c-d) show the projected horizontal and vertical emittances, and (e-f) shows both the kinetic energy (left 
axis) and rms bunch length (right axis). 



allowed users to explore adjustments of the injector op- near real time with measurements. Using this interface, 
tics in simulation while in the control room, often in and starting from a single optimized setting of the injec- 
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TABLE II. 






Element 


Parameter 


19 


pC/bunch Values 


77 pC/bunch Values 


laser 


pinhole (mm), intensity cut off 


(%) 


1, 40 


2, 35 


laser 


rms pulse length (ps) 




8 


8 


DC gun 


voltage (kV) 




350 


350 


solenoid 1 


peak field (T) 




0.032 


0.031 


buncher 


voltage (kV), phase (deg) 




50, -90 


60, -90 


solenoid 2 


peak field (T) 




-0.020 


-0.020 


SRF cavity 1 


voltage (kV), phase (deg) 




1491, -10 


1491, -10 


SRF cavity 2 


voltage (kV), phase (deg) 




1953, -16 


1953, -7 


SRF cavity 3 


voltage (kV), phase (deg) 




1386, 


1386, 


SRF cavity 4 


voltage (kV), phase (deg) 




1386, 


1386, 


SRF cavity 5 


voltage (kV), phase (deg) 




1386, 


1500, -20 


A3 quad 1 


field gradient integral ([T/m]- 


m) 


0.013 


0.013 


A3 quad 2 


field gradient integral ([T/m]- 


m) 


-0.033 


-0.033 


A3 quad 3 


field gradient integral ([T/m]- 


m) 


-0.016 


-0.016 


A3 quad 4 


field gradient integral ([T/m]- 


m) 


0.029 


0.029 


Bl quads 


field gradient integral ([T/m]- 


m) 


-0.017 


-0.016 



TABLE III. Measured and simulated projected horizontal (a), and vertical (b) emittances. 

(a) Horizontal projected emittance data. 



19 pC Measurement Type 


en,4lOO%) 


en,490%) 


en,x(core) 


/core 


en,a; (core) //core 


Projected EMS 


0.33 ±0.02 /im 


0.23 ±0.02 /im 


0.14 ±0.01 /xm 


67% 


0.21 ±0.01 /xm 


Time-resolved EMS 


0.28 ±0.02 /im 


0.21 ±0.01 /xm 


0.14 ±0.01 /xm 


72% 


0.19 ±0.01 /xm 


GPT Simulation 


0.31 /im 


0.19 /im 


0.07± /xm 


59% 


0.12 /xm 


77 pC Measurement Type 


en,a. (100%) 


en,490%) 


en,x{coie) 


/core 


Cn, a: (core) //core 


Projected EMS 


0.69 ±0.05 fim 


0.51 ± 0.04 /im 


0.28 ±0.2 /xm 


64% 


0.44 ± 0.03 /xm 


Time-resolved EMS 


0.66 ±0.05 fim 


0.48 ± 0.04 /xm 


0.29 ± 0.2 /xm 


67% 


0.43 ± 0.03 /xm 


GPT Simulation 


0.72 /xm 


0.44 /xm 


0.17 /xm 


51% 


0.33 /xm 





(b) 


Vertical projected emittance data. 






19 pC Measurement Type 


en,y (100%) 


6n,^(90%) 


en, y (core) 


/core 


en,y(cOre)//core 


Projected EMS 
GPT Simulation 


0.20 ±0.01 /xm 
0.16 /xm 


0.14 ±0.01 /xm 
0.11 /xm 


0.09 ±0.01 /xm 
0.06 /xm 


70% 
64% 


0.13 ±0.01 /xm 
0.09 /xm 


77 pC Measurement Type 


en,y (100%) 


en,y{90%) 


en, y (core) 


/core 


en,y(cOre)//core 


Projected EMS 
GPT Simulation 


0.40 ± 0.03 /xm 
0.37 /xm 


0.29 ± 0.02 /xm 
0.25 /xm 


0.19 ±0.01 /xm 
0.11 /xm 


70% 
59% 


0.27 ±0.01 /xm 
0.19 /xm 



tor model, optics sets for both 19 pC and 77 pC bunch 
charges were found which kept both the simulated rms 
beam sizes small throughout the injector, in addition to 
preserving the minimized emittance at the end of the 
merger section. These settings were then loaded into the 
injector, and the phase space data for each bunch charge 
was taken. 

The resulting data sets include the vertical and hor- 
izontal projected phase spaces, as well as the time- 
resolved horizontal phase space at the merger EMS, and 
the energy spread distribution in the A4 straight sec- 



tion. Overall, we found excellent agreement between 
measurement and simulation. For both bunch charges, 
the agreement between the measured projected 100 and 
90% emittance values was within 30% of the simulated 
values in both transverse planes. We point out that for 
77 pC/bunch, the measured 90% emittance in vertical 
plane, as well as the core emittance in both planes, meets 
the ERL design specification of e^ < 0.3 /xm. The pro- 
jected emittance in both transverse planes demonstrates 
the correct scaling with bunch charge shown in Eq. (IT]). 
Using the time-resolved horizontal phase data, the longi- 



17 




40 60 

f (%) 

(a) Horizontal 19 pC/bunch data. 




(b) Horizontal 77 pC/bunch data. 
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(c) Vertical 19 pC/bunch data. 
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(d) Vertical 77 pC/bunch data. 



FIG. 14. Emittance versus fraction curves: (a-b) the curves computed from the measured horizontal phase space data at 19 
and 77 pC/bunch, (c-d) the curves computed from the measured vertical phase space data at 19 and 77 pC/bunch. 





TABLE IV. Simulated and measured rms energy spread as a function of bunch charge. 


Data Type 


19 pC/bunch, A4 section 77 pC/bunch, A4 section 


19 pC/bunch, Bl merger 77 pC/bunch, Bl merger 


GPT Simulation 
Measured 


0.16% 0.27% 
0.14 ±0.01% 0.26 ±0.01% 


0.12% 0.21% 

NA NA 



tudinal bunch profile and time-resolved emittance were 
computed. The measured rms bunch length for both 
bunch charges was at or below the 3 ps specification, and 
agreed with simulation to within 5%. For both bunch 
charges, the time-resolved core emittance met the ERL 
specification. Finally, an estimation of the energy spread 
of the beam in merger was found by measuring the energy 
spread in the straight section. Agreement between the 



measured and simulated rms energy spread was within 
13% for both bunch charges. 

These results represent a significant advancement in 
high-brightness photoinjectors. The measured emit- 
tances in this work set a new record low for DC photoin- 
jectors producing beams with comparable bunch charge. 
To put these results in a broader picture, it is instruc- 
tive to compare the performance of the Cornell injector 
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OpC 



19 pC 



77 pC 




data 



GPT 



1 1.5 

X (mm) 

(a) Horizontal phase-space as a function of bunch charge. 

OpC 19 pC 77 pC 




data 



GPT 



y (mm) 
(b) Vertical phase-space as a function of bunch charge. 



FIG. 15. Comparison of the measured and simulated projected transverse phase-space as a function of bunch charge. Plot (a) 
shows the the horizontal phase-space, while (b) shows the vertical phase-space. Corresponding emittance values can be found 
in Table Uni 



for its designed application of a 5 GeV x-ray ERL to the 
beam quality of existing storage rings. For this compari- 
son, we assume a 100 mA, 1 nm-rad horizontal emittance 
storage ring with 10~^ energy spread and 1% coupling 



factor, representing the best of existing third generation 
light sources \^3^. As a figure of merit for non- Gaussian 
beams, it is convenient to use the effective transverse av- 
erage beam brightness over the rms energy spread of the 
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GPT: a^ = 2.2 ps, data: o^ = 2.1 ps 



GPT: a^ = 3.1 ps, data: o^ = 3.0 ps 





Time (ps) 

(a) Instantaneous current at 19 pC/bunch. 
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(b) Instantaneous current at 77 pC/bunch. 



FIG. 16. Comparison of the measured beam current to GPT simulation. The estimated uncertainty in the rms bunch lengths 
was (a) ±0.1 ps and (b) ±0.2 ps. 
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(a) 19 pC/bunch slice emittance data. 
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(b) 77 pC/bunch slice emittance data. 



FIG. 17. The 90% and core slice emittance for: (a) 19 pC/bunch, (b) 77 pC/bunch. 



beam at the location of an undulator: 



Jx ' Jy 



^x\Jx) ' ^y\Jy) 



crs 



(8) 



Here Cx and Cy are the transverse geometric emittance 
values as a function of the horizontal and vertical beam 
fractions, respectively. The energy spread is included in 
this expression to reflect the fact that undulators with 
larger number of periods can be more efficiently utilized 
for beams with smaller energy spread. In an ERL, the 
energy spread after the main linac will be defined by the 
RF curvature and the bunch length according to (27r/rf • 
G-tf' I\f2 [44 . Using our 19 pC/bunch data, and assuming 
the full repetition rate, the estimated energy spread and 



effective average brightness of a 1.3 GHz, 5 GeV ERL 
yields a higher transverse brightness over the best storage 
ring by a factor of 20. 

Looking forward, we point out that the measurements 
shown here demonstrate two crucial points: (i) that low 
emittances reported previously in simulations [71 18] are 
well within the reach of the next planned iteration of the 
photoinjector; and (ii) the relevant physics and control 
parameters required to produce these low emittances are 
now understood. In developing a plan for reducing the 
emittance further, we note that optimization results in- 
dicate that lower emittances and shorter bunch lengths 
at the end of the merger are possible at higher beam en- 
ergies [30]. As a result, the optimal photoinjector for a 
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(a) 60% isosurface rendering of the 19 pC time-resolved phase-space, (b) 60% isosurface rendering of 19 pC sHce emittance data showing 

inner isosurfaces. 



X (mm) n. . p 

-0.05 omO£L^^ 



yySxXlOOO 




X (mm) „ 1 „ 

-0.05 0.00 0^J^_5i2 



ry^xxiooo 



Time (ps) 



(c) 60% Isosurface plot of the 77 pC time-resolved phase-space. 
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(d) 60% isosurface rendering of 77 pC slice emittance data 
showing inner isosurfaces. 



FIG. 18. Three dimensional plotting of the time resolved phase-space 



future ERL light source will operate at higher beam en- 
ergies (roughly 12 MeV) than those used in this work [2]. 
Eq. (fTl) shows two more directions for further improve- 
ment. For a given bunch charge, the emittance in this 
equation can be reduced by lowering photocathode MTE, 
or by increasing the accelerating field at the cathode. In 
fact, the results for the vertical emittance demonstrate 
that the emittance in this plane is dominated by the ther- 
mal emittance, and thus colder cathodes are required. 
Currently, there is an active cathode research program at 
Cornell University dedicated to improving cathode per- 
formance [45]. Already, cathodes with MTE values as 



low as 30 meV have been experimentally realized both 
for negative affinity and multi-alkali photocathodes [46] . 
In parallel, Cornell is developing an improved DC gun, 
in order to overcome the current voltage limitation. The 
new gun design features a segmented insulator with guard 
rings [47] in order to minimize damaging the insulator 
from field emission. Lastly, improved laser shaping will 
aid in creating bunches with more linear space charge 
fields. According to the rough scaling law in Eq. (IT]), 
as well as more detailed calculations reported in [2J [8] , 
these improvements are expected to reduce the emittance 
in the photoinjector by roughly a factor of 3, resulting in 
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a beam brightness roughly 10 times higher than reported 
here. This abihty to independently improve critical ele- 
ments in the injector, resulting in better performance, is 
one of the major strengths of a linac based accelerator. 

We acknowledge Mike Billing for his help with non- 
linear BPM corrections, as well as Bill Schaff for provid- 
ing the GaAs photocathode used in these measurements. 
This work is supported by National Science Foundation 
(Grant No. DMR-0807731). 



Emittance Definitions 

Here we provide the relevant emittance definitions used 
in this work to describe non-Gaussian phase space distri- 
butions. We use the standard definition of the normalized 
transverse rms emittance: 
1 



mc 



V{x'^){pI) - {xp^y 



VWWFW^¥^iW. 



(9) 



where 7 and Px are the normalized energy and trans- 
verse velocity of each electron. In this and all fol- 
lowing expression the subscript "n" is used to distin- 
guished between the normalized emittance e^ and geo- 
metric emittance e, which are related by e^ = (7/^) • e. 
In this and all subsequent expressions, {u) denotes the 
average over the particle distribution in phase space: 
{u) = JJ u{x^Px)p{x^Px)dxdpx^ where p{x^Px) is the nor- 
malized 2D phase space distribution function. The rms 
emittance as a function of beam fraction is defined as fol- 
lows [48 . For an area in phase space 7ra, an ellipse with 
Twiss parameters given by 




is defined so that the phase space region enclosed by the 
ellipse is given by d{a^T) = {x : x^T~^x < a}, where 
X = (x^Px)^ . The twiss parameters in T are varied until 
the fraction of particles enclosed in the ellipse is maxi- 
mized. Labeling this phase space region D{a)^ the beam 
fraction is defined as: 



/(a) = max { j I p{x,Px)dxdpx 

ld{a,T) 



p{x,px)dxdpx. 



(10) 



D{a) 



The corresponding fractional emittance takes the form 

^n{0) = :^\/{x^)d{pI)d - {xPx?D. (11) 

where {u)d = j^ JJ^^^^up{x,Px)dxdpx. The para- 
metric curve defined by {/(a), €n{a)} is the emittance 
vs. fraction curve Cnif)- Also important for understand- 
ing emittances of non- Gaussian beams are the definitions 



of the core emittance and corresponding core fraction 

[lasg: 



en (core) 



dcn 
df 



•) /core • ^nv/corej — ^nvCOrCj. 



/^o 



(12) 



For comparison purposes, the emittance vs. fraction 
curve for 2D uniform, elliptical, and Gaussian distribu- 
tions have been computed. To do so the correlation be- 
tween X and Px has been removed and the coordinates 
rescaled so that the distributions can be written as radial 
functions of the normalized coordinate f = \/^^^^^- 
Additionally, the distributions are parameterized so that 
the resulting emittance vs. fraction curve en(/) is nor- 
malized: en{f = 1) = 1. Fig 19 (a) I shows each of the 
three distributions as a function of the normalized radial 
coordinate. The corresponding emittance vs. fraction 
curves are shown in Fig. |19(b)[ From these curves the 
90% and core emittance (relative to the 100% emittance) 
can be computed. Table \V\ gives these ratios, as well as 
the core fraction, for each distribution. For additional 
discussion on the connection between core emittance and 
brightness, see |48] . 
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TABLE V. For comparison the corresponding scaled emit- 
tance and fraction for various phase space distributions are 
given. 



Distribution Type 



£n(9Q%) Cnjc 



U 



£n(core)//c( 



uniform 
elHptical 
Gaussian 



0.90 
0.87 
0.74 



1 100% 

5/6 87% 
1/2 72% 



1 
0.96 
0.69 



